A bosonic Josephson junction controlled by a single trapped ion 
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We theoretically investigate the properties of a double- well bosonic Josephson junction coupled to 
a single trapped ion. We find that the coupling between the wells can be controlled by the internal 
state of the ion, which can be used for studying mesoscopic entanglement between the two systems 
and to measure their interaction with high precision. As a particular example we consider a single 
87 Rb atom and a small Bose-Einstein condensate controlled by a single 171 Yb + ion. We calculate 
inter-well coupling rates reaching hundreds of Hz, while the state dependence amounts to tens of Hz 
for plausible values of the currently unknown s-wave scattering length between the atom and the 
ion. The analysis shows that it is possible to induce either the self-trapping or the tunneling regime, 
depending on the internal state of the ion. This enables the generation of large scale ion-atomic 
wavepacket entanglement within current technology. 

PACS numbers: 03.75.Gg, 03.75.Lm, 37.10.Ty, 34.50.Cx 



The phenomenon of tunneling in a Josephson junction 
is one of the most striking quantum effects that can be 
observed on a macroscopic scale. It has been studied in 
superconductors 

0,11 , superfluid helium 0,13] 

and ultra- 
cold atomic gases 5 8]. In the latter, the interaction 
between the particles plays a crucial role, leading to phe- 
nomena such as (quantum) self-trapping A plethora 
of studies have been performed on this system con- 
sidering, for instance, one-dimensional junctions d, [To| 
and analogs of photon assisted tunneling [ll], EH ■ The 
bosonic Josephson junction may be well suited to study 
macroscopic entanglement (l3l. Il4j. like in Schrddinger- 
cat-type superpositions that are instrumental in quan- 
tum information science and in the study of the crossover 
from quantum to classical physics. Additionally, highly 
non-classical states are an important resource to measure 
external (classical) fields with high precision [l5| . 

Recently, systems involving a combination of ultra- 
cold quantum gases and trapped ions have become avail- 
able (m [I?} . These experiments have great potential to 
study atom- ion collisions in the quantum regime pH , [l9j , 
cold chemistry (2(J [2l[ , and to explore sym pat hetic cool- 
ing of ions by means of clouds of cold atoms [22( . Further- 
more, accurately controlled trapped ions could be used 
to perform in situ measurements of cold atomic gases in 
optical lattices [23| , or to induce a ^m-sized bubble in a 
bosonic Tonks-Girardeau gas [24| . 

In this work we investigate a double well bosonic 
Josephson junction that is controlled by a trapped single 
ion as schematically shown in Fig. 1. The electric field 
of the ion polarizes the atoms leading to a long range at- 
tractive potential. For small distances between the atoms 
and the ion, however, the interaction is strongly repulsive 
and depends on the relative spin orientation of the atom 
and ion via the singlet and triplet scattering length (25j . 



Hence, the tunneling is controlled by the trapped ion, 
especially by its internal (qubit) states {| f), | 4-)}. This 
is of great interest, as it enables the generation of meso- 
scopic entanglement between atomic systems and ions 
and to study their interaction with high precision. 

Our strategy to attack the hybrid Josephson junction 
problem consists of two steps: first, we solve the problem 
at the single-particle level, where only a single atom in- 
teracting with the ion is trapped in the double well poten- 
tial. We will show that there is a strict relation between 
the inter-well coupling and the state-dependent atom-ion 
interaction. Moreover, we show that this dependence is 
rather large such that it could be observed experimen- 
tally. Then, the single-particle solutions will be used 
to study the many-body problem within the two-mode 
approximation. Finally, we shall analyze the impact of 
inelastic collisions and heating, and discuss experimental 
implementations . 

We start with the single-particle problem. The first 
ingredient is the interaction between an atom and an ion 
caused by an induced atomic dipole due to the electric 
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FIG. 1. (Color online) Bose-Einstein condensate trapped in 
a double well potential with an ion in the center. The internal 
state of the ion can be tuned by laser light and controles the 
tunneling rate J. 
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FIG. 2. (Color online) a) Wave functions for the ground states 
$ g (z) (blue) and $ e (z) (red), at q = 2.06iT; the black curve 
shows the potential with u q = 2tt 3.9 kHz. The wavefunctions 
are multiplied by a constant for clarity, b) Spectrum for a 
double well as a function of q and short-range phases of <f> = 
— 7r/4 (solid lines) and cj> = — n/3 (dashed lines) with the 
ground states shown in bold red and blue. The dashed vertical 
line marks the separation q = 2.06-R*. 



S 150 : 
x 

b 100 



5 50 

0- 



0.4 0.5 
0&r) 



FIG. 3. Inter- well coupling J for various short-range phases 
with ujg = 2tt 3.9 kHz and q = 2.06R*. 



field of the ion. At large distances it is given by: 
lim V ia (r) = 



(1) 



with C4 = e 2 a p /2, where e is the charge of the ion, 
and a p is the static polarizability of the atom. The 
potential (|T|) is characterised by the length scale R* = 
v /2 j uC 4 /ft 2 and the energy scale E* = K 2 / (2/.i(R*) 2 ), 
with [i the reduced mass. For r — ^ 0, Eq. ([T]) does not 
hold anymore and the potential becomes strongly repul- 
sive. The exact form of the potential in this regime is gen- 
erally not known. Given this, we will use quantum defect 
theory, where all the information on the short-range in- 
teraction is put in a single quantum defect parameter, the 
so-called short-range phase <p, which is related to the ion- 
atom s-wave scattering length aj a = — R* cot(</>) [25L |2(| . 

Now we consider the experimentally realistic case in 
which the ion is trapped much more tightly than the 
atom, cjj 3> uJai with u)^ a the trapping frequency of the 
ion and atom. We also assume that the ion motion is 
cooled to the ground state (e.g., by sideband cooling). 
These assumptions enable us to make an approximation 
in which the ion is treated statically at the position = 
0. We further assume that the ion and atom cannot 
undergo spin changing collisions, allowing the use of a 
single channel model. This is justified when the singlet 
and triplet scattering lengths between the ion and atom 
are similar [26j . or when there is only one open channel 
for the considered states. Under these conditions the 
Schrodinger equation for the atomic wave function, in 
units of R* and E* , is given by: 



Ei/)(t) 



^V 2 



+ V dw (v)--j)^(r), (2) 



where Vd w ( r ) is the atomic double well potential, and m a 
the atomic mass. We consider a double well of the form: 



V dw (r) 



E*q l 



(z 2 -q 2 ) 2 + a^( X 2 + y 2 ), (3) 



which has minima at r 



ping frequency uj q = 



(0, 0, ±q) with local trap- 
and b is the barrier height 



between the wells. For convenience we assume that 
Lu q = uj ai with uj a the atomic trapping frequency in the 
(x, ?/)-directions. We consider varying the inter-well sepa- 
ration q without changing the local trapping parameters: 
only the inter-well barrier height is changed. Besides this, 
a = (R*/l ) 4 and l = ^/h/m a uj a . 

For r — > the atomic trapping potential and en- 
ergy are negligible compared to the atom-ion interaction 
and the solution can be obtained analytically: ipi{r) = 
N^/r [Ji+i/ 2 (0 +tan(<S)^+ 1/2 (£)]. We use this solu- 
tion as a boundary condition for numerically solving 
Eq. ([2) [27[. Here, Ji+i/ 2 (0 and ^ + i/ 2 (0 arc spherical 
Bessel functions with orbital angular momentum quan- 
tum number I, £ = ^Jm a / ' \ir and N is a normalization 
constant. The mixing angle S is related to the short- 
range phase: 5 = —(f) — Itt/2. In our study we neglect 
small corrections due to the /-dependence of <f>, because 
in the ultracold collisional regime basically only s-wave 
scattering occurs [2|| . 

As a specific example to our model, we consider a single 
trapped 87 Rb atom in the \F a — 2,m,F a = 2) state and 
a 171 Yb + ion with logical states | 4-) = \F% = 0, iriFi = 0) 
and I 'f) = \Fi = l^mpi = 1)- When the ion is in state 
I f), the dynamics is described by the triplet scattering 
length at- When the ion is in the state | J.) the interaction 
depends on both the singlet and triplet scattering length. 
Since these values are not currently known, we rather 
assume plausible values for the corresponding short-range 
phases <fi^ and <f>±. 

We set the parameter a to 10, corresponding to u> q = 
2ir 3.9 kHz. For the chosen atom and ion i?*=0.306 ^m 
and E* jh = 2n 935 Hz. In Fig. [2] the spectrum as a 
function of the inter-well separation q is shown. For 
large q, the spectrum is equivalent to that of a har- 
monic oscillator, since there is no coupling between the 
wells. For smaller q the levels split in two as the inter- 
well coupling increases. The ground state splits in the 
symmetric (ground) and anti-symmetric (excited) super- 
positions $ 9;e (r) = ($i(r) ± &n(r))/y/2 with energies 
Eg and E e as shown in Fig. [2] Here, $fl,i(r) denotes 
the left/right localized wavepackets. When the splitting 
AE = E e — Eg is much smaller than the energy gap to 
the other levels in the spectrum, a two-mode approxima- 
tion can be employed, with the inter- well coupling given 
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FIG. 4. (Color online) a) Time dependence of the inter-well 
distance q for generating state dependent tunneling, b) Prob- 
ability distributions at various times. The initial state pre- 
pared at q — 2.3K* is an equal superposition of the two ground 
states ip(z,t = 0) = ($ g (z) — <& e («))/\/2. At the end of the 
sequence the atom has tunneled when the ion is prepared in 
| t) (top panel b) and returned to the left well after tunneling 
multiple times when the ion is prepared in | 4-) (lowest panel) . 
For the calculation we assume <j>f — —n/4 and <j>i — — n/3. 
At the end of the sequence > 98% is in the right/left well 
respectively. 



calculation of this process is shown in Fig. |4Jd), where 
the potential and the atomic probability distributions are 
plotted. The setup can also be thought of as an interfer- 
ometer, in which the entangling-disentangling dynamics 
can be monitored to study the interaction between the 
ion and the atom. Similar studies have been performed 
e.g. with cold Rb atoms in optical lattices [28j . where 
differences in scattering length for different spin states 
lead to such dynamics. 

Now let us consider the many-particle scenario. To 
this aim we apply a two-mode approximation in which 
the quantum field operator is expanded in the local 
modes §>l,r(t) of the double well potential, that is, 
V>(r,t) = § L (r)c L + $ R (r)c R , where Cj (cf) with i = L,R 
is the annihilation (creation) operator of an atom in the 
left or right well. The local modes are the single-particle 
ground states of the system that we have obtained al- 
ready. The two-mode approximation provides a reason- 
able description of the dynamics when many-body in- 
teractions yield only slight modifications to the ground 
state. This requires N <C lo/a a , with N the number 
of atoms, Iq = (Ixlylz) 1 ^ 3 the size of the ground state 
wave packet, and a a the atomic s-wave scattering length. 
Under these assumptions, the system can be effectively 
described by the following Bosc-Hubbard (BH) Hamilto- 
nian 0,111: 



by J = AE/(2h). As Fig. [2|d) shows, this coupling de- 
pends on the internal state of the ion. In this example 
we chose 4>i = — tt/3 and <fr^ = —tt/4. At the position 
q = 2.0QR*, the coupling amounts to Jj- = 2ir 11.2 Hz 
and Ji = 2ir 39.2 Hz. Without an ion the tunneling 
rate would be only 2ir 0.96 Hz. In Fig. |3l the inter-well 
coupling for various other short-range phases </> is shown. 
This analysis clearly shows that there is a relatively large 
dependence of the inter- well coupling on the short-range 
phase, which should be experimentally observable, and 
J depends on the qubit state via </>. Hence, we can treat 
the inter-well coupling as an operator and write it as: 
J = ^4.1 I) (I I + Jf\ t)(t I- Moreover, typical spectra 
show much larger energy gaps to neighbouring states, 
such that the two-mode approximation is well justified 
in a many-body calculation. 

We can use the state dependence in the double-well 
to design an experimental sequence that entangles the 
internal state of the ion with the position of the atom. 
To this aim, we introduce time dependence into the well 
separation q(t) [Fig. H]a)] and solve the time-dependent 
Schrodingcr equation numerically for the sequence: 1) 
The ion is prepared in the internal state ( | 4-) + 1 1})/ V% 
and the atom is initially trapped in the left well. Then 
we adiabatically bring the wells together. 2) The atom 
tunnels and we wait until the atom is either in the right 
or left well, depending on the (internal) state of the ion. 
3) We take apart the wells again such that tunneling 
stops. 4) The resulting atom-ion state is the entangled 
state \ip final) oc | 1)\$r) + e if} \ f)\$ L ), where the phase 
j3 depends on the details of the sequence. A detailed 



Hbh = hja + hUfi , 



(4) 



where a 

at - 



c r cl is the tunneling operator, and 

n = (c r cr — c j l cl)/2 describes the population imbalance 
between the two wells. Because of the spin-dependence of 
the atom-ion interaction we also write the on-site atom- 
atom interaction as an operator U. 

We recall that the above outlined two-mode model en- 
ables to identify, depending on the ratio U/J, three dif- 
ferent regimes Q: the"Rabi" regime (i.e., weak atom- 
atom interaction) ; the intermediate " Josephson" regime, 
and the "Fock" regime (i.e., strong atom-atom interac- 
tion). In our study we are interested in the Rabi and 
Fock regimes. We show that depending on the internal 
state of the ion, it is possible to induce either "Rabi" 
oscillations between the two wells or to keep the con- 
densate trapped in a single well (i.e., the self-trapping 
regime). It is useful to introduce the interaction param- 
eter A = UN/2 J. Self-trapping regime is predicted to 
occur for A > A c = 2 [H- 

Now the on-site interaction energy can be estimated by 
using the single-particle wave functions leading to U — 
E*!f- j dr|$ L (r)| 4 with U = 8Tra af i/m a R* [2g. Here 
we neglected terms such as E*^ J dr|$L(r)| 2 |(f> iR (r)| 2 , 
which arc much smaller than U. Then, in order to quan- 
tify the quality of the description in localized wavepack- 
ets, which is crucial for the BH model, we calculate the 
probability e of finding the atoms in the 'wrong' region, 
namely e = J R |<j>i(r)| 2 dr 3 , and the integral taken over 
R, the right half space, (z > 0). 

Let us consider a weaker trap, such that both the Rabi 
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FIG. 5. (Color online) Relative population in the left well 
as a function of time for | f) (blue) and | \) (red) N = 20 
atoms. Obtained by integrating Eq. Q numerically, assuming 
(j>f = — 7r/4 and tf>± — ir/4. 



and the self-trapping regimes are achieved for a rela- 
tively large number of atoms. By choosing ui a = uo q = 
2?r 200 Hz (a = 0.026) and q = 6.39 R*, the two-mode 
approximation requires N <C 100 and we find the follow- 
ing parameters: Uf = 2tt 0.9 Hz, J-f- = 2tt 1.7 Hz for 
cfo = — 7r/4 (a ia = R*) with e = 0.004. Whereas for neg- 
ative scattering lengths a; a = —R*, i.e., 4>i = 7r/4, we 
have Ui = 2tt 1.0 Hz, Jj. = 2tt 42.7 Hz with e = 0.017. 
For N — 20 atoms we have Aj = 5.2 (i.e., self-trapping 
regime) and A;t = 0.23 (i.e., the Rabi regime). The re- 
sulting BH dynamics is shown in Fig. [5J which clearly 
shows that large scale entanglement can be achieved. For 
N = 100 atoms, we would find Aj = 25.8 and A^: = 1.17. 
For <p t = tt/3, i.e. (a ia = -R*/y/3), we have for N = 20 
atoms A| = 2 and for N = 100 atoms Aj = 9.8. 

We underscore, however, that in the self-trapping 
regime the BH-model is in disagreement with the full 
numerical many-body calculation especially at long time 
scales [30] . The time scales are set by the Rabi period 
tRabi = 71*/ J. With the above outlined parameters we 

find: 4,m/*b-w = 25 ( 8 ) fOT <h = "^4 (tt/3). There- 
fore, we might expect reasonable agreement between the 

Bose-Hubbard model and a full quantum many body cal- 

i~ i 
culation at short times, that is, t -C t Rabi but t ~ ^ a hi 

and not too far into the self-trapping regime, as in Fig. [5] 
We will now discuss some of the limitations of our 
model and possible experimental issues. As we already 
noted before, the atom and ion could change their spin 
states in a collision (2|| . When the ion is in the state | f) , 
spin changing collisions cannot happen due to the con- 
servation of the quantum number Mp = m a + m, = 3: 
there are no other states with Mp = 3. When the ion is 
in the state | \), Mp = 2, we have to consider the mix- 
ing with the states \F.i,mpi, F a ,m Fa ) = 1 1111) , 1 1121) 
and 1 1022). However, all these states have significantly 
higher dissociation energies, since the hypcrfine splitting 
of m Yb+ (12.6 GHz) is much larger than that of 87 Rb 
(6.8 GHz). Thus, only resonances with bound states can 
occur, which have a much lower density than trap states. 
Consequently, we can reasonably assume that also for 
|0022) no spin flips occur. 

Concerning inelastic collisions, it is important to note 
that our proposal requires low atomic densities. We es- 
timate the peak density of the atomic wavepackcts to be 
about 10 20 /m 3 . Now, assuming that many-body interac- 



tions do not perturb significantly the form of the wave- 
function, this would lead to loss rates of a few Hz fl6j . 
Secondly, in our proposal the atoms and ion are not di- 
rectly overlapped, the ion rather perturbs the atomic 
wavefunction while it is held in its own trapping poten- 
tial, reducing the density and terefore the inelastic scat- 
tering a factor >10 more. This setup may also reduce 
heating due to the 'micromotion' of the ion in its time- 
dependent trap [3l], HH- However, even in the presence 
of moderate inelastic collision rates, our proposal could 
be implemented, since these effects are easily separated 
from experimental errors: one would simply observe that 
there is no (Yb + ) ion at the end of the experiment. 

Regarding cooling and heating mechanisms, we first 
note that the ion can be cooled to the ground state 
with effectively unit efficiency using sideband coolings. 
However, given the time scales required for the proposed 
experiments, heating may cause problems. In state of 
the art ion traps, such as in [33l |. heating rates of a few 
phonons/s are achieved, such that an experiment on the 
tens of ms timescale should be possible. The coherence 
of the atomic wavepacket can be sustained for tens of ms 
even for well-separations of several /im [|| . 

Finally we mention that Josephson junctions for cold 
atoms have been implemented in optical dipole traps and 
in micro- fabricated magnetic traps Such systems 

can be combined with ion trap technology. One interest- 
ing approach would be to integrate atom and ion traps 
in a single micro-trap. Recent experiments have demon- 
strated micro fabricated wires or permanent magnets on 
ion surface traps [3~4T - f37| , with the purpose of performing 
quantum logic gates without lasers by magnetic coupling 
to the ions. In Refs. [3^, magnetic field gradients of 
2-20 T/m are reported which would make it possible to 
magnetically trap 87 Rb atoms with trapping frequencies 
reaching 20 kHz, 50 ^m above the trap surface. 

In conclusion, we studied the tunneling dynamics of 
a single atom and a small BEC in combination with a 
trapped ion. We found that the junction transmission 
can be controlled by the internal state of the ion mak- 
ing it possible to engineer large ion-matterwave entan- 
gled states. The system could be scaled up by consider- 
ing Josephson junctions coupled by strings of ions. The 
ion-atom hybrid quantum system may lead to precision 
measurement of external fields, or collisional properties 
where short-range interactions may be measured by de- 
termining tunneling frequencies. 
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I. ADDITIONAL MATERIAL 
A. Diagonalisation of the double well hamiltonian. 

To numerically solve Eq. (2) in the paper, we make use 
of renormalized Numerov method \3S with 



as a boundary condition. Our approach is to obtain a 
set of basis functions if>°(r) with energies E° in which 
to expand the solutions to the double-well problem. The 
superscript is put to mean that there is no inter-well 
separation, as will become clear below. A convenient 
choice for this purpose is an isotropic atomic harmonic 
oscillator V(r) = ar 2 /j/m„ with a single ion trapped in 
its center, where a = (R*/Iq) 4 and la = y/ K/ m a oj a is 
the size of the atomic ground state. The solutions ip° 
have the form tp° = Y^ipni (r) / r, where Y™ are spherical 
harmonics and ip n i (r) / r denotes the radial wavefunctions 
to be found. We obtain a set of 1250 cigenfunctions and 
energies including bound as well as trap states [26[. The 
procedure is similar to the one used in Ref. (26J. 
Now we consider the double well of the form: 



V dw (r) 



E*q 4 

2 + ! 



a — (as + y ) 

m„. 



o// 



The first term in Eq. ([6]) corresponds to the harmonic po- 
tential used for determining the basis functions ^(r), 
and therefore this part of the Hamiltonian is diago- 
nal. On the other hand, for the remaining terms in 
Eq. (|6]), "0nz( r ) arc not cigenfunctions, and thereby we 
have to compute the corresponding matrix elements (us- 
ing z — r cos 9). We use the expansion onto the solutions 
of the isotropic harmonic well t/>(r) = J^k ^V'fcWi where 
k denotes the pair of quantum numbers (n, I). Since the 
potential does not depend on the azimuthal direction, 



the quantum number m is conserved and, without loss 
of generality, we set m = 0. Both the wavefunctions and 
the corresponding energies are obtained by diagonalising 
the Hamiltonian, whose matrix elements are given by: 



Hhk> = 



El 



fxaq 



3pia (2 
Okk' - Ml 



2m a 



kk' 



4m a q 2 kk ' ' 



Mr) = Ny/r [J|+i/ a (0 + taa(*)^ +1/2 (f)] (5) M ^ 



dr s tpp(r) cos J (6»)r J ' V°(r), 



(6) 



with j = 2, 4. The clear advantage of this procedure is 
that the matrix elements do not depend on q, such that 
spectra are easily obtained as a function of the inter- well 
separation. 

B. Short range state dependence 

For r — > the interaction between atom and ion is 
described by the quantum numbers \FMpIS) [25|, with 
F = F a + F, , the total nuclear spin I = i a + U and to- 
tal electron spin S = s a + Sj. Depending on whether 
S = or S = 1, the interaction is described by the sin- 
glet or triplet scattering length a s and at, respectively. 
The quantum states for large r, labeled by the hyper- 
fine states \Fi,m,Fi, F a ,mFa) are connected to the quan- 
tum numbers for the short-range interaction via a frame 
transformation 25|. When the ion is in state | f), the 
dynamics is described by the triplet scattering length a t , 
since S — 1. On the other hand, when the ion is in the 
state | 4-) , the dynamics depends both on the singlet and 
triplet scattering length, since this state is a linear com- 
bination of short-range states \FMpIS), in particular we 
have: 




|F a = 2 ! m Fo =2)U) = -i|2220)+y| 



which can be obtained by construction of the frame trans- 
formation in terms of Clebsch-Gordan coefficients and 
Wigner 9-j symbols (see Ref. (HI). Here, the right hand 
side is labeled according to \FMpIS). 
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